Entangled Bessel beams 
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Abstract: Orbital angular momentum (OAM) entanglement is investi- 
gated in the Bessel-Gauss (BG) basis. Having a readily adjustable radial 
scale, BG modes provide a more favourable basis for OAM entanglement 
over Laguerre-Gaussian (LG) modes. The OAM bandwidth in terms of BG 
modes can be increased by selection of particular radial modes and leads 
to a flattening of the spectrum. The flattening of the spectrum allows for 
higher entanglement. We demonstrate increased entanglement in terms of 
BG modes by performing a Bell-type experiment and violating the appro- 
priate Clauser Home Shimony Holt (CHSH) inequality. In addition, we 
reconstruct the quantum state of BG modes entangled in high-dimensions. 
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1. Introduction 

The use of the orbital angular momentum (OAM) eigenstates of photons in quantum informa- 
tion science became attractive after it was shown, theoretically QEl and experimentally 0, 
that OAM is conserved during spontaneous parametric down-conversion (SPDC) — the az- 
imuthal indices of a pair of down-converted photons add up to that of the pump beam. As a 
result, the pair of down-converted photons are naturally entangled in terms of their azimuthal 
indices or OAM. Entanglement is a desirable property in quantum information applications, 
where it has been used in quantum ghost imaging [4], quantum cryptography [6l and in 
quantum computing algorithms Q. Any transverse spatial modal basis defines an infinite- 
dimensional Hilbert space, allowing more information per photon to be used in quantum in- 
formation applications. 

Although the Laguerre-Gaussian (LG) modes are currently the popular choice for theoretical 
analyses of quantum information systems based on OAM, the modal basis used in quantum 
information experiments are seldom pure LG modes; see Ref. 0. The optical manipulation 
and detection of LG modes requires control of the radial index, which governs the radial shape 
of each LG mode's intensity profile. We shall refer to the modes that are detected without this 
radial control as azimuthal modes. The azimuthal index is not associated with a unique quantum 
state, but with a subspace of the total Hilbert space. If the radial part of the mode changes, a 
photon propagating through an optical system may suffer a loss of its quantum entanglement, 
even if the azimuthal index is unaffected. As the complexity of quantum information systems 
increases, the requirements for the fidelity of the quantum states can be expected to become 
more demanding. Although the connection between the spatial modal profile of an optical beam 
and OAM was initially made with specific reference to LG beams [9], the same property applies 
to any optical beam with a rotationally symmetric intensity profile. As a result, Bessel beams 
GoHH) and Bessel-Gaussian (BG) beams llT2l also have quantized amounts of OAM associated 
with them. Each photon in such a beam also carries an amount of OAM equal to lh, where £ is 
the azimuthal index of the mode. 

In this paper, we investigate the use of BG modes fT2l as a basis for photonic quantum 
information experiments. The radial index associated with LG modes is replaced by a contin- 
uous scaling parameter for the radial part of the BG modes. The relative ease and accuracy 
with which one can manipulate the azimuthal and radial dependencies of BG modes make 
them a favourable choice for high-fidelity quantum information systems. We compare the high- 
dimensional density matrices measured in the BG basis to those measured in the azimuthal 
basis. Our results indicate that under the correct experimental conditions, the generation and 
measurement of high-dimensional entanglement in the BG basis is advantageous compared to 
the azimuthal basis. Our results indicate a possible advantage in using BG modes over azimuthal 
modes for both generating and measuring OAM entanglement. 



2. Bessel-Gauss Modes 



The electric field of a BG mode is given by 
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where £ is the azimuthal (mode) index (a signed integer); J^(-) is the Bessel function of the 
first kind; k r and k z are the transverse and longitudinal wave numbers respectively, obeying 
k 2 + k 2 = k 2 . The initial radius of the Gaussian profile is COo and the Rayleigh range is zr = 
7TC0q I X . It is clear from Eq. (Q]) that the radial component of the mode can be scaled by altering 

k r . 

The quantum state for the photon pairs produced in SPDC can be written in terms of BG 
modes, 

W = L// MkriM) \£Ms\-ZM)idk rl dk r2 , (2) 

where \d£(k r \,k r 2)\ 2 denotes the probability amplitude for measuring a signal photon in state 
\£,k r \) s and an idler photon in state | — £,k r 2)i. The radial scaling parameters k r \ and k r 2 in 
Eq. © can be adjusted in a continuous manner to optimize the bandwidth of the OAM spec- 
trum. 

For particular BG modes specified by k r , the two-photon state for d dimensions can be written 

as 

W=IXM>i. (3) 

t 

where £ ranges over d different values. The sets of ^-values for different dimensions d are 
chosen as in Ref . lfT3l . 

The density matrix of a two-photon quJit quantum state can be reconstructed in terms of a 
Bloch representation, 

d 2 -i 

P = b m,n?m®?ni (4) 
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where b mA are complex coefficients with Z?o,o = l/d 2 for normalisation; T p are the d- 
dimensional generalised Gell-Mann matrices for p = 1 ... (d 2 — 1 ) ; To is the d-dimensional 
identity matrix. These coefficients are determined by performing a tomographically complete 
set of measurements. An example of an over-complete set consists of the BG states \£,k r ) to- 
gether with the superpositions of two BG states 

|a) = -^[|^i,t r ) + exp(ia)|^,*r)]. (5) 

Here, a is the phase between the modes of the superposition states and k r is chosen as a fixed 
value for all basis states. The quantum state in Eq. © consists of two quJits, which results in a 
d 2 x d 2 density matrix and thus requires at least d 4 — 1 measurements to match the number of 
matrix elements lH4l . However, if an over-complete set of measurements is chosen, additional 
information is available to compensate for measurement errors due to natural photon number 
fluctuations. An uncertainty of 5% was observed for the measured coincidence counts. 



3. Experimental Setup 

Our experimental setup is shown in Fig. [T] A mode-locked laser source with a wavelength 
of 355 nm and an average power of 350 mW was used to pump a 0.5-mm-thick type I BBO 
crystal to produce non-collinear, degenerate entangled photon pairs via SPDC. The plane of the 
crystal was imaged onto two separate spatial light modulators (SLMs). The SLMs were used 
to select the particular pair of BG modes that were to be detected, as explained below. The 
SLM planes were re-imaged and coupled into single-mode fibers to extract the pure Gaussian 
modal components. The fibers were connected to avalanche photodiodes, which then registered 
the photon pairs via a coincidence counter. The single count rates and were recorded 
simultaneously with the coincidence count rates Q and were accumulated over a 10 second 
integration time. 
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Fig. 1. (a) Experimental setup used to detect the OAM eigenstate after SPDC. The plane of 
the crystal was imaged onto two separate SLMs using lenses, fi = 400 mm and f2 = 750 
mm, where the BG modes were selected. Lenses f$ = 300 mm and £4 = 1.5 mm were used 
to image the SLM planes through 10 nm bandwidth interference filters (IF) to the inputs 
of the single-mode fibres. Examples of a phase-only binary bessel hologram with a helical 
phase of £ = 2 for different values of k r are shown in (b) k r = rad/mm, (c) k r = 2\ rad/mm 
and (d) k r = 35 rad/mm. The inset shows a back-projected CCD image of a binary bessel 
mode with helical phase £ = 2 and radial wavevector k r = 2l rad/mm measured at the plane 
of the crystal. 

In order to select BG modes using SLMs, we encode a binary Bessel phase function to 
approximate the Bessel function, thereby avoiding the need to perform complex amplitude 
modulation (15), which results in the loss of optical power. The binary Bessel function is then 
combined with a helical phase function to determine the azimuthal index t of the detected BG 



mode. The phase-only hologram used to project the photons into BG modes is given by the 
transmission function 

T(r) = sign{J^ r r)}exp(^0), (6) 

where sign{-} denotes the sign-function. Spatial filtering introduced in each of the detection 
arms can be used to remove the binary nature of the transmission function, causing the ef- 
fective modes that are detected to be the BG modes. Without the binary Bessel function, the 
transmission function is the same helical phase function that is traditionally used to measure the 
azimuthal modes. However, each OAM mode with a particular value of \t) is a superposition of 
different radial modes that have the same azimuthal index. The binary Bessel function, which is 
encoded on the SLMs, enables us to select one particular radial mode from this superposition. 
As a result of selecting one radial mode, we anticipate two notable features in the measured 
coincidence counts. Firstly, the spiral bandwidth of the BG spectrum will change depending 
on the chosen radial mode. Secondly, as only one radial mode from the superposition is se- 
lected, the coincidence count rate will be lower when the Bessel function sign{/^} is included 
compared to when it is not. The shape of the spectrum will also change as a result of the bi- 
nary Bessel function, because one only measures the part (one-dimensional slice) of the full 
two-dimensional spectrum that corresponds to a particular value of k r . 

Following the method derived in Ref. fT6lh we used superpositions of OAM states with ±£, 
in analogy to the polarisation states, to demonstrate the violation of a Bell-type inequality. The 
superposition states for BG modes given by Eq. (0) were oriented at angles and a# for the 
signal and idler photons, respectively. 

Using the method described in Ref. fT3lh we calculate the density matrix for higher dimen- 
sions by minimizing the Chi- square quantity, 
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where p are the normalized coincident counts, which are interpreted as experimentally meas- 

(p) 

ured probabilities, and p\ ' are the predicted probabilities calculated from the guessed density 
matrix. We construct the density matrix so that it is both Hermitian and positive semi-definite 
with unitary trace lfT4l . This is necessary as the effects of statistical variations of the coin- 
cidence count rates and experimental inaccuracies commonly produce reconstructed density 
matrices with negative eigenvalues and a non-unitary trace. 



4. Experimental Results 

In the first experiment we vary the radial wave number k r on each SLM and monitor the coin- 
cidence counts for each £ mode. For k r = the phase functions on the SLMs consist of only the 
helical phase. As a result the k r = measurements correspond to the conventional azimuthal 
modal spectrum. The OAM bandwidth can be quantified as the standard deviation (square root 
of the second moment) in I, given by 
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or as the Schmidt number ifTTIl 
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Both quantities are plotted as a function of k r in Fig.|2](a). 



(a) (b) 




Fig. 2. (a) OAM bandwidth (blue) and Schmidt number (red) are shown as a function of 
k r . As the radial wavevector increases from k r = rad/mm (azimuthal modes), the OAM 
bandwidth increases, (b) OAM spectrum for azimuthal modes (red) and BG modes for 
k r = 21 rad/mm (purple) and k r = 35 rad/mm (green), in terms of the modal weightings. 
The maximum coincidence count rate is lower for both BG modes, however a broader 
spectrum is observed with increasing k r . 



We observe almost a doubling in the OAM bandwidth from k r = rad/mm (azimuthal modes) 
to k r = 35 rad/mm (BG modes). This tuneability in the radial component of the field allows for 
the OAM spectrum to be customised through application of Eq. $2&. In Fig.[2](b) one can ob- 
serve the broadening of the OAM spectrum as the k r is increased from k r = rad/mm (red) 
to k r = 21 rad/mm (purple) to k r = 35 rad/mm (green). The maximum coincidence count rate 
decreases for BG modes, but its spectrum is broader than that of the azimuthal modes. The 
full- width half-maximum (FWHM) widths of the spectra indicate that the BG modal spectrum 
at k r = 35 rad/mm is almost twice as wide as the azimuthal modal spectrum obtained at k r = 0. 
The wider bandwidth can be explained as an effect of the more selective measurement of the 
transverse spatial modes, thanks to the specification of the radial definition of the mode (Eq.0). 
As a result one observes a reduction in the central region of the spectrum. Note that the reduc- 
tion in the central region of the spectrum does not correspond to an attenuation over all radial 
modes, but rather to a more selective detection of a particular radial mode. The reduction in 
the central region of the spectrum makes Procrustean filtering in higher dimensional entangle- 
ment experiments fT8l superfluous. Beyond a certain t- value the BG spectrum coincides with 
the azimuthal modal spectrum. This can be understood as a consequence of having a large £- 
value with a fixed width Gaussian envelope: the Gaussian suppression of the function beyond 
a certain radius causes the definition of the radial dependence of the mode for large I- values to 
have little or no effect. Although the spatial resolution of the SLMs would become a limiting 
factor, larger values of k r are expected to give OAM spectra with larger bandwidths. The wider 
bandwidth allows more basis states to be incorporated in quantum information systems. 

To verify that the higher-dimensional quantum state produced in the SPDC process is en- 
tangled, we use it to test a Bell-type inequality. Clauser, Home, Shimony and Holt (CHSH) 
provided a method to experimentally test Bell's theorem, stating that no local hidden- variable 
theory could simulate the kind of statistical predictions made by quantum mechanics fT9l . We 
fixed oca from Eq. © at four different angles: a A = 0°,45°,90°, 135°, allowing a# to vary from 
0° to 180°. The Bell parameter, as defined in Ref. |[T6lL is given by 



S = E(a A ,a B ) - E(a A , a' B ) +E(a' Al a B ) - E(a! A , a! B ), 



(10) 



Angle of hologram B, a E [degrees] 



Fig. 3. Sinusoidal behaviour of the normalised coincidence counts as a function of the 
angular position of the holograms, for £ = 1 at positions a A = 0° (blue), 45° (pink), 90° 
(green) and 135° (yellow). The insets show the holograms used for oc A = 0° and oc A = 45°, 
where the phase varies from (black) to K (white). 



where 



E(a A ,a B ) 



C(a A , a B ) + C(a A + f , a B + f ) - C(a A + %,a B )- C(a A ,a B + ^) 
C(a A , a 5 ) + C(a A + ^ , a# + + C(a A + a B ) + C(a A , a 5 + ^) ' 



(11) 



with C(a A , a B ) being the coincidence count rate for the particular orientation of each hologram. 
We see a clear violation of the CHSH inequality for BG modes (k r = 21 rad/mm) — our Bell 
curves are shown in Fig. [3] For the two-dimensional BG subspace of £ = 1, we find the Bell 
parameter in the CHSH inequality to be S = 2.78 ± 0.05 > 2, which violates the inequality by 
15 standard deviations. 

We also reconstructed high-dimensional two-photon density matrices in the BG basis for 
k r = 21 rad/mm using a full quantum state tomography ll20l . The density matrices for d = 2 
and d = 5 are shown in Fig. (4). The purity of the state Tr(p 2 ) was calculated as 0.981, where 
its closeness to 1 indicates that the observed quantum state is close to a pure state. To quantify 
the mixture of the measured state we calculate the linear entropy I2H Sl = 4/3 [1 — Tr(p 2 )] . We 
find the linear entropy of our quantum state to be Sl2 = 0.06, indicating that our state is close 
to a pure state (the linear entropy for a pure state is zero, whereas a completely mixed state is 
charaterised by a linear entropy of one). 

We also calculate the fidelity 



which is a measure of how close our reconstructed state is to the target state pr- In our case 
the target state is the (pure) maximally entangled state. We find a fidelity of F = 0.96. These 
results are indicative of entanglement of our BG modes. We performed high-dimensional state 
tomography measurements from d = 2 to d = 8 and use the result to compute the fidelity and 
the linear entropy, as indications of entanglement and purity, respectively. We found that the 
higher-dimensional quantum states in terms of the BG modes for a particular k r have higher 
fidelity than the quantum states in terms of azimuthal modes, while remaining relatively pure, 
as shown in Fig. Uc) and (d). These measurements are compared with the threshold states, 
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Fig. 4. Results from a full quantum state tomoography of a BG mode with k r = 21 rad/mm. 
(a) & (b) Graphically representation of the real part of the density matrix for dimension, 
d = 2 and d = 5, respectively, (c) Linear entropy and (d) fidelity as a function of dimension. 
The red triangles represent the measured data for the azimuthal modes, the green squares 
represent the measured data for the BG modes and the blue circles represent the threshold 
states in Eq. JT3l >. 



which are given by 

which lie on the threshold of the high-dimensional Bell inequality l22l . Here p™ m is the prob- 
ability above which the Bell inequality is violated and is the identity matrix of dimension d 2 . 
The fidelity values ranged from F2 = 0.96 to F$ = 0.79 for the BG modes, while the values for 
linear entropy were found to be Sl2 = 0.06 and Sus = 0.35. 

5. Discussion and Conclusion 

We have successfully demonstrated that BG modes are entangled in the high-dimensional OAM 
degree of freedom. Crucially, the ability to select the radial component allows us to record BG 
states with a very broad spiral spectrum. The resulting high-dimensionally entangled states 
are closer to the appropriate maximally entangled state compared to the case where no radial 
selection is performed. Such control of the experimental conditions results in a higher informa- 
tion capacity per photon pair and allows for more states to be utilized in quantum information 
processes. The only cost of the increased spiral spectrum of the recorded modes is a lower 
coincidence count rate, which is due to the required modal selectivity. Due to the particular 
interference properties of BG modes, entangled modes in this basis may offer increased ro- 
bustness when propagating through turbulent media. The finite space-bandwidth product of the 
SLMs limits the achievable spectrum of the BG modes; however, alternative methods may offer 
further improvements in OAM entanglement with BG modes. 



